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Abstract 

In this paper, we investigate the probabilistic characteristics of a unit driven by 
Dichotomous Markov Noise (DMN), as an external random life increasing and decreas- 
ing shocks. Using DMN, we will define two new aging classes of the overall increas- 
ing/decreasing (OIL/ODL) nature in the long time behavior, which are separated by 
an exponential steady state regime. In addition, a moment inequality is derived for 
the system whose life distribution is in an overall life decreasing (ODL) class. We use 
this inequality to devise a nonparametric testing procedure for exponentiality against 
an alternative overall decreasing life distribution. 

General Theory The statistical theory of positive aging systems is an interdisciplinary 
domain of research, which is encountered in many physical, chemical, biological, and engi- 
neering safety systems. Here, we study the positive aging process driven by a Dichotomous 
Markov Noise (DMN). Describing the deterioration or positive aging of life in the engi- 
neering and biological systems, based on different aging criteria, has been the subject of 
investigation for several decades. For more technical details we refer the reader to Mar- 
shall & Olkin (2010), Shaked & Shanthikumar (2010), Rolski (1975), Barlow & Proschan 
(1981), Deshpande et al. (1986), Abouammoh (1988), Mitra & Ansi (2008), and Ibaranim 
& Sepehrifar (2009). The plan of this paper is as follows: In section one, we introduce 
the concept of Dichotomous Markov Noise (DMN) and introduce two new aging classes of 
Overall Decreasing Life (ODL) and Overall Increasing Life (OIL). The moment inequality 
approach is introduced in section two. In section three, a new statistics test, based on 
U-statistics theory is established to test exponentiality against an alternative distribution 
F, where F belongs to an overall decreasing life (ODL) class of distributions. In section 
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four, Monte Carlo simulations for sample sizes n = 10, 20, 30 are conducted for this test. 
The calculated empirical power of the test shows the excellent power of this test for some 
common alternative distributions. The novelty in this study is the development of a non- 
parametric testing procedure for testing expoenetiality against the two new aging classes 
for a unit driven by Dichotomous Markov Noise (DMN). Furthermore, this proposed test 
is of special interest to biostatisticians, who are typically contend and are limited to data 
that comes at successive patient visits to a facility; such data ignores crucial events that 
may occur between visits. 



1 Introduction 

Let X{t), a continuous positive random variable, describing the lifetime (aging) of a mechan- 
ical unit, living organism or some other physical, electrical or biological system (Deshpande 
et al. , 1986). We interpret aging as the time that the system has worked satisfactorily 
without experiencing a failure and then enters into a failed state (Deshpande et al. , 1986; 
Kochar & Wiens, 1987). In this context, a system with the older age has a shorter lifetime 
(Kochar & Wiens, 1987). In addition, we assume that X{t) is described by Dichotomous 
Markov stochastic process. 



1.1 Dichotomous Markov Noise (DMN) 

First order systems driven by dichotomous noise are of practical interest in a wide variety of 
biological, physical and engineering problems (Bena, 2006; Bicout, 1997; Yarahmadian et al. 
, 2011). Similar to an Ornstein-Uhlenbeck process, DMN is colored noise and defined as a 
two- valued stochastic process with state space values v± and constant transition frequencies 
of A±; the increase (+) and decrease rate (— ) for \± occur with probabilities p±{t). The 
switching of v{t) are Poisson processes. The evolution of the unit age, X{t), is represented 
by the following system of differential equations (Bena, 2006; Yarahmadian et al. , 2011): 



(^l ^ (p+\ _ (-^+ ^- \ IP+ 



dt \P-J V ^+ ~^-J \P 

Using the initial conditions p+{0) = 1 and p-{0) = 0, we can represent the solutions as 
follows: 



where p+{t) +p-{t) = 1 for all times, and j = ^^^x represents the mean time between 
switches of v{t), otherwise known as the rate relaxation time. 
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Figure 1: Dichotomous Markov Noise. 

Remark 1.1. DMN is a time-homogeneous Markov process and is completely characterized 
by the following transition probabilities: 

(1.4) = Pr{vit) = vMO) = Vj), hi e {±} 

The temporal evolution of transition probabilities are written as: 

.^ ± (P-jit)\ _ (-K \- \ (P-j{t)\ 

with the following solutions: 

(P—{t) P^+{t)\ _ 1 (\+ + A_e-(^++^-)* A+(l - e-(^++^-)*) 
^'■^^ \P+-{t) P++{t)) A+ + A_ VA-(l-e-(^++^-)*) A_+A+e-(^++^-)* 

In the stationary case we have: 
(1.7) Pr{v = v+)= ^ , Pr{v = v^) = ^ 
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Lemma 1.2. The aging life of the system is characterized by V = ^'^x++x ■ 
Proof. We define the instantaneous average rate as: 

(1.8) v(t)=v^p+{t)-v-p-{t) 
Using (1.2) and (1.2), we get 

(1.9) v{t) = V{1 - e-(^++^")*) + ^;+e-(^++^-)* 
and 



(1.10) V = Urn v{t) 



□ 

Definition 1.3. When V < 0, life decreasing shocks are more frequent and the system is 
in an Overall Decreasing Life (ODL) class, whereas when V > life increasing shocks are 
more frequent and the system is in an Overall Increasing Life ( OIL ) class. 

The following lemma captures the situation when the life distribution of the system falls 
in the steady-state: 

Lemma 1.4. In steady-state, the mean unit age and its probability distribution function 
(pdf), are obtained through an exponential distribution: 

(1.11) A= Z^'^" , m = \e-'^ 

Proof. The time evolution of the unit age distribution P±{x,t) is represented by: 

n ION ± (P+{^,t)\ _ f-v+9^-X+ A_ \ fP+{x,t) 

^ ' dt\P-{^,t)J V A+ y_d_-x.j \P.{x,t) 

In the steady state, since there is no dependence on time , the time derivatives are set 
to zero and (1.12) is reduced to the following: 

n 1o^ ± fP+{x,t)\ _ (-^ ^\ (P+{x,t) 

^ ^ dx[P-{x,t))-[-^ ^)[P-{x,t) 

or equivalently to: 

Using the fact that lima;^oo -P±(a;) = 0, the solutions to (1.13) are written as: 
(1.15) P±{x) = A±e~i 

where A± are two real constants. Now f{x) = P+{x) + P-{x) and f{x)dx = 1 will 
establish the result. □ 
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2 Statistical Approach: Moment Inequality Approach 



Let F{x) be the c.d.f. of X{t) and let F{x) = 1 - F{x) = Pr[X{t) > x] be the survival 
function. The survival function of a unit of age xq, Fxq{x), is defined as the conditional 
probability that a unit of age xq (i.e., a unit having the initial age xq) will survive for an 
additional x units of time: 



(2.1) F{x\xo) = Pr[X{t)>x + xo\X{t)>xo]-- + 



F{xo) 



If Xq = 0, then by definition F{x\0) = F{x) is the survival function of a new unit. Any 
study concerning the phenomenon of aging must be based on F{x\xq) and functions related 
to it. Associated with X is the notion of " random remaining life" at age xq, denoted 
by XxQ which has a survival function ^p^^^^ .The concept of mean remaining life (MRL) 
has found applications in biometry, actuarial science, and engineering, nonparametric life 
testing and reliability analysis. For reviews of the applications of MRL, the reader can refer 
to Muth (Muth, 1977), Hall and Wellner (Hah & Wellner, 1981), Gupta and Gupta (Gupta 
&: Gupta, 1983) and the various references cited in these works. 

Now, let Xi,X2, ■■■ be a sequence of independent identically distributed random lives with 
survival function F = 1 — F and finite mean A ( A is described by (1.11)). 

Remark 2.1. The distribution (1.11) is meaningful if and only if = . 

Remark 2.2. Abouammoh et. al. (Abouammoh et al. , 1993) introduced a model of accu- 
mulated additive damages that cause the failure of a device when it exceeds a certain value 
threshold. Abouammoh et. al. (Abouammoh al. , 1994) ^^■^o established an empirical test 
statistics for testing exponentiality when the alternative to exponentiality belongs to the new 
better than renewal used in expectation (NBRUE). 

Let N[t) be the total number of shocks governed by Poisson process, which take place 
in [0, t). Then the remaining life at t is L{t) = ^^^^^^ Xi — t, assuming J2i^=i^ Xi — x > 
where x is the fixed threshold. In this setting L{t) converges weekly to a renewal random 
life denoted by X with survival function W{x) = ^ F{u)du (Ibrahim &: Mugdadi, 2010). 
Note that X can be studied through the remaining (residual) life X^q = [X — xq\X > xq\ 
at age xq. 

Theorem 2.3. For all x > 0, W{x) = \E{X - x)I{X - x) where I{u) = 1, if u > and 
otherwise. 

Proof see (Ibrahim & Mugdadi, 2010). □ 

Definition 2.4. A life distribution F on (0,oo), with F{0~) = is called overall decreasing 
life (ODL), if 

/oo 
W{x)dx < AW{t), t > 0, 



5 



where A = F{u)du < oo. 



Remark 2.5. Life distribution F belongs to ODL if and only if < j^. 

There are a number of classes that have been suggested in the hterature to categorize 
distributions based on their aging properties or their dual. Among the most practical 
aging classes of life distributions, we study the Increasing Failure Rate (IFR) . A survival 
distribution is said to be in the IFR class if and only if ^^p^°'^ is decreasing in x for all 
xq. The definition says that the probability of an individual of age x survives an additional 
xq period of time is decreasing with time. 

Remark 2.6. Barlow and Proschan (Barlow & Proschan, 1981) showed that an IFR life 
distribution F with mean A is bounded below by e~ a for t < K. 

Lemma 2.7. The implications between the IFR and ODL classes are as follows: 

IFR =^ ODL 

Proof IFR ^ ODL: 

This follows from direct application of Remark 3.4 and Remark 3.5. 

□ 

Theorem 2.8. // a lifetime distribution F belongs to the ODL class, then for all integers 

r > 

(2.3) 

E\—^—X^+^Xi — Xr+^X2+ , ^ < AE\—XiMlf ^—XiAir+A 

l2(r + l) 1 2 r + 2 ^ 2(r + 3) ^ J" lr + 1 ^'^ r + 2 ^-M 

where Xi, X2 are two stochastically independent copies of life X and Mi^2 = min{Xi,X2). 

Proof. F belongs to the ODL class if W{x)dx < AW{t) . Apply Theorem 2.3, one may 
rewrite (2.2) as follows: 

(2.4) / / t'"F{t)W{x)dxdt < A / fF{t)W{t)dt 
Jo Jt Jo 

L.H.S= / / fF{t)W{x)dxdt 
Jo Jt 

(2.5) =iMr(H-'"'^^+H*} 
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r-oo 

R.H.S = A / fF{t)W{t)dt 
Jo 

Umin(Xi,X2) _ 
{Xi - t)I{Xi > t)fF{t)dtj 

= e\^—Xim:+^ — —Xim;+^\ 

□ 

When r = we deduce the following special case: 
Corollary 2.9. Setting r = in (2.3); 

(2.7) e[^X^xI - \xlX^ + \x\\ < AE{XiM,,2 - \mI^} 

3 Statistical Testing of X Aging (ODL Alternative) 

The goal of this section is to create a symmetric kernel, allowing the application of U- 
statistics theory for testing exponentiality against any other alternatives in the ODL aging 
class. We use corollary (2.9) to test i^o : is exponential vs Ha : F is ODL. We define the 
measure of departure from the null hypothesis as: 

(3.1) S = E[^X^Xi - \xlX2 + \xl - A(XiMi,2 - \m12)\- 
To make the test statistic scale invariate, we adjust 5 by and define 

(3.2) A = ^e\^-X^xI - \xlX^ + \xl - A(XiMi,2 - ^Mf ^)} 
Note that A = under the null hypothesis and A < under the alternative. 

Let X\^X2^ . . .Xn denote a random sample from F. A U-statistics form of b is defined 

as 

(3-3) ^=^^7;^ E ^^^i^-\^i^3^\n-nM,,^\x,Mi^^ 

i<i<j<n 
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and 



(2.6) 



Now, we can estimate A by A = -p-, where X = ^ Yli -^i- To perform this test, we calculate 
A from data and reject Hq in favor of Hi if the normal variate value Zi-a exceeds a^^^na A. 
In this calculation one may find the value of variance through the standard argument of 
[/-statistics theory (Lee, 1990) as follows: 



Theorem 3.1. As n ^ oo, y/n{A — Aq) is asymptotically distributed normal with mean 
and variance 



-^Var{-X! + x^dFix) - -XfF{Xi) + -{X^ x^dF{x) 

+ / x^dF{x) - Xl / xdF{x)) \ 
Jo Jo ^ 



(3.4) ^ 



under the null hypothesis (Tq = 1.173. 
Proof. Setting 

(3.5) HXi,X2) = {IxiXl - \xlX2 + \xl - XlMi^2 + \xiMl^] 
and then applying the standard theorem of U-statistics (Lee, 1990), we obtain 

(3.6) a2 = yar{E{$(Xi,X2|Xi)} + ^{$(X2,Xi|Xi)}}, 

in which 
(3.7) 

£;{Xi2min(Xi,X2)|Xi} = ^jXi^ min(Xi, X2)I(Xi < X2)\Xi^ + E^XlX2l{Xi > X2)|Xi} 

= XlF{Xi) +Xf [ ' xdF{x) 
Jo 

and 
(3.8) 

^{x|min(X2,Xi)|Xi} = ^{x|min(X2,Xi)/(Xi < X2)|Xi} + min(X2, Xi)/(Xi > X2)|Xi} 



/•oo 

Xl / {x^ + x'^)dF{x) 
Jxi 



Direct calculations will establish the results. □ 

4 Simulated Power Studies 

The empirical power estimates of our test statistics are obtained based on 10000 replicates 
of the Monte Carlo method simulating common alternative life distributions with increasing 
failure rate. Weibull, Gamma, and linear failure rate family with the distribution functions, 

(4.1) Fi{x) = e-^'", X > 0, 9>1, 
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Table 1: Power Estimates Using a = 0.05 



Distribution 


Parameter 9 Sample Size 

77 — in 77— 20 77 — Sn 


Gamma Family 


1 0.4693 0.7078 0.7870 

2 0.8635 0.9794 0.9914 

3 0.9641 0.9980 0.9993 


Linear Failure Rate Family 


1 0.7672 0.9700 0.9920 

2 0.8455 0.9890 0.9977 

3 0.8767 0.9929 0.9999 


Weibull Family 


1 0.4641 0.7093 0.7754 

2 0.9930 ^1 

3 ssl ^1 ^1 



(4.2) F2(a;) = e-^-^', x > 0, ^ > 0, 
and 

(4.3) F3(x) = — — / e""n^~^du, x > 0, e>0 
respectively, are among these alternatives. 

We report the power values for samples of size n = 10, n = 20, and n = 30 in Table 1 when 
1 < 6. These estimates of power were obtained for tests where a = 0.05. More extensive 
simulations and results, along with the tables, are available upon request from the authors. 
We note from Table 1 that is 9 >1 (increasing the distance of the exponential), the proposed 
test get better power; on the other hand the power decreases as 6 approaches to 1 , which 
reduces to the ordinary exponential distribution. 

4.1 Numerical Example 

Abouammoh et.al (Abouammoh et al. , 1994) worked on data which represented 40 patients 
suffering from leukemia from a Ministry of Health Hospital in Saudi Arabia. The life time 
number of days, in increasing order, are: 115, 181, 255, 418, 441, 461, 516, 739, 743, 789, 
807, 865, 924, 983, 1024, 1062,1063, 1165, 1191, 1222, 1222, 1251, 1277, 1290, 1357, 1369, 
1408, 1455, 1478, 1549, 1578, 1578,1599, 1603, 1605, 1696, 1735, 1799, 1815, 1852. 
The test statistic for this data set is A = —0.871605. When a = 0.05, (TQ^^n'ih. = 
—3.615245 <^ —1.6448. Thus, from this calculated value, we must reject the null hypothesis 
and we conclude that this specific set of data has the ODL property. 

5 Conclusion 

The main goal of this paper is to apply nonparametric testing procedures to test if the life 
time of a system, driven by DMN, falls in the overall decreasing life (ODL) class or not. The 
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mathematical tools and the proposed test introduced here are of general application and 
should be of general interest to those individuals studying biostatistics, biology, and safety 
management systems engineering. The results of this paper can also be easily extended to 
situations involving right-censored data. It is our hope that this method and its application 
will open the way to the study of more complicated and realistic examples. 
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